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ABSTRACT 



Context. The dominating reionization source in the young universe has not yet been identified. Possible candidates include metal poor dwarf 
galaxies with starburst properties. 

Aims. We selected an extreme starburst dwarf, the Blue Compact Galaxy Haro 11, with the aim of determining the Lyman continuum escape 
fraction from UV spectroscopy. 

■ Methods. Spectra of Haro 1 1 were obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE). A weak signal shortwards of the Lyman 
' break is identified as Lyman continuum (LyC) emission escaping from the ongoing starburst. From profile fitting to weak metal lines we 
derive column densities of the low ionization species. Adopting a metallicity typical of the H II regions of Haro 11, these data correspond to a 
hydrogen column density of ~10 19 cirT 2 . This relatively high value indicates that most of the LyC photons escape through transparent holes in 
the interstellar medium. We then use spectral evolutionary models to constrain the escape fraction of the produced LyC photons. 
Results. Assuming a normal Salpeter initial mass function we obtain a Lyman continuum escape fraction of f esc ~ 4-10%. We argue that in a 
hierarchical galaxy formation scenario, the upper limit we derive for the escape rate allows for a substantial contribution to cosmic reionization 
by starburst dwarf galaxies at high redshifts. 
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1. Introduction 

According to WMAP data, the epoch of reionization started at 
a redshift of z ~ 20 (Kogut et al. I2003i . Observations of the 
high opacity in the Lyman continuum (hereafter LyC) emis- 
sion of galaxies (Becker et al., 120011 Fan et al. |?002 1 and 
the results from the recent studies of the intergalactic medium 
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(IGM) around SDSS quasars (e.g. Mesinger & Haiman 2004 1 
indicate that the final stage of the epoch of reionization may 
have occurred at redshifts z ~6-7. There is still no consensus 
about what could have been the major ionization sources of 
the IGM during the reionization epoch - stars, normal galax- 
ies, AGNs or some exotic mechanism (Yan & Windhorst 2004 
Meiksin & White 2003 1. The rapid decline of the space den- 
sity of quasars at redshifts z >3 (Madau et al. 119991 Fan et al. 
2001 1 makes it likely that their contribution to the LyC produc- 
tion rate required for the reionization is insignificant, although 
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it has been proposed that miniquasars fed by intermediate-mass 
black holes may have been important (Madau et al. 2004). 

At lower redshifts, observations of a sample of Lyman- 
break galaxies at z ~ 3.4 (Steidel et al. 2001 1 seem to show that 
galaxies could account for a significant fraction of the photon 
flux required for reionization. The observations indicate a con- 
tribution of LyC photons a factor of ~5 higher than for quasars 
at the same epoch. These results have however not been con- 
firmed by other similar studies. On the contrary, it has been 
argued (Fernandez-Soto et al. 2003) that a maximum of 4% 
of the LyC flux could escape from luminous galaxies at red- 
shifts 1.9< z <3.5. This is also close to the upper limit of the 
escape fraction of low redshift galaxies (Leitherer et al. 1995 
Hurwitz et al. 119971 Deharveng et al. 1200 H Heckman et al. 
2001 1. The lack of a strong trend with redshift indicates that 
galaxies at higher redshifts also have low escape fractions, in- 
sufficient to keep the intergalactic gas ionized. Alternative ion- 
ization sources that have been discussed are massive pregalac- 
tic Pop III stars and young globular clusters. Most probably the 
reionization is ruled by different sources at different redshifts. It 
is therefore important to establish observational constraints on 
the LyC fluxes from the various candidates at different epochs. 

Our information about the galaxy luminosity/mass function 
at the epoch of reionization is very meager, and it is not at all 
clear what mass scale would dominate the LyC photon produc- 
tion rate. In a galaxy formation scenario where mergers are of 
major importance one would expect starburst dwarf galaxies to 
be frequent among the first generation of galaxies. Indeed the 
number density of so called Luminous Compact Blue Galaxies 
(LCBGs) appears to increase with redshift (Lilly et al. 1998 
Mallen-Ornelas et al. U999> . It is therefore important to inves- 
tigate if this type of galaxy could have been an important ion- 
ization source at high redshift. Dwarf galaxies are however too 
faint to be observed at high redshifts. We are therefore depen- 
dent on theoretical models and observations of nearby galaxies 
that can be used as templates of the first generation of dwarf 
galaxies. Recent models do not give much hope for Lyman 
escape from galaxies of masses larger than 10 7 solar masses 
(Kitayama et al. 2004) but these predictions have to be tested 
observationally. 

We have searched for local dwarfs with properties similar 
to young dwarf galaxies, that could be used as test objects. A 
few such studies have already been carried out but only pro- 
vided upper limits to the escape fraction (e.g., Giallongo et al. 
[T9971 Devriendt et al. fT99l Shull et al. fT999l . However, pre- 
dicted equivalent widths of the Balmer emission lines in star- 
bursts are in some cases larger than observed (e.g. Bergvall 
[T9851 Breso lin et al.. [T9991 Mas-H esse & Kunth [T9991 Moy 
et al. 120011 Stasinska et al. 2001). Several possible explana- 
tions to this were discussed by Mas-Hesse & Kunth C1999). 
One possibility is that the most massive stars remain hidden in 
their parental clouds until most of the fuel has been consumed. 
This would mimic a truncation of the IMF at high masses. 
Mas-Hesse and Kunth also discussed the possibility that a large 
fraction of the ionizing photons are absorbed by dust in the H 
II regions. To obtain agreement with the observations, a large 
amount of the photons, ~ 60% would have to be absorbed. This 
possibility has recently been investigated by Hirashita et al. 



(2003 ) where the escape fraction of the LyC photons from the 
H II region itself is estimated from model comparisons. Based 
on input data from a sample of star forming galaxies they find 
that, in the mean about 40% of the photons are trapped, which 
is less than what would be needed according to Mas-Hesse and 
Kunth. But the scatter is high so in some cases this certainly 
would be a possibility. Hirashita et al. seem to find that the 
escape fraction is higher is galaxies with properties similar to 
Haro 1 1 (e.g. the BCG ESO 338-IG04). Therefore an alterna- 
tive, or complementary, explanation for the low flux of the re- 
combination lines, may be that a significant fraction of the LyC 
photons is in fact escaping from the galaxy. This is the possi- 
bility that we study in this paper. 

H II regions in nearby galaxies appear to leak 40-50% of 
the LyC photons produced by the hot stars into the diffuse 
interstellar medium (Ferguson et al. 1996 Oey et al. 1997). 
Although it cannot be completely excluded that the gas may 
be ionized by single O stars located in situ, a comparison 
of the observed emission line strengths with theoretical mod- 
els (Iglesias-Paramo & Munoz-Tunon 2002 , Wood & Mathis 
2004 ) give support to the leakage interpretation. An even more 
extreme case of a density bounded H II region around a mas- 
sive star cluster was reported by Leitherer et al. ( 1996). They 
estimate that about 75% of the LyC radiation is leaking into 
the diffuse ambient medium. Depending on the distribution of 
the gas, the location where these photons are absorbed may be 
very different from case to case. In some situations the diffuse 
ionized gas may be much more extended than the H II region 
producing the photons and will be ignored by the observer who 
consequently will underestimate the star formation rate. 

The other alternative for the discrepancy between theory 
and observations of the Balmer emission line strengths is that 
the dust is unevenly distributed across the starburst region so 
that the younger regions are more reddened than the older ones 
(Calzetti et al. I1994> . This could reduce the equivalent width 
of Ho- with up to a factor of « 2. The only way to really dis- 
tinguish between a global leakage into to intergalactic medium 
and the other two alternatives is to search for leaking photons 
shortward of the Lyman limit. 

In previous studies of starburst dwarfs (Bergvall 1985 
Bergvall & Ostlin 120021 Ostlin et al. IT9991 l200ll we have 
focused on a few luminous (-19>Mg >-21) metal poor blue 
compact galaxies (BCGs). They have properties similar to the 
LBGs at intermediate redshifts (e.g. Guzman et al. 2003). Here 
we report on such a study using FUSE, the Far Ultraviolet 
Spectroscopic Explorer. 

1.1. The target galaxy 

The target galaxy of our choice is Haro 1 1 (=ESO 350-IG38) at 
ff20oo=00 /, 36 m 52.5*, <5 2 ooo=-33 49'49". From the emission-line 
spectrum we derive a heliocentric radial velocity of Vq=6180 
km s _1 , corresponding to a distance of 86 Mpc, assuming a 
Hubble constant of 72 km s _1 Mpc -1 . We have previously stud- 
ied this galaxy over a wide frequency range, from UV to radio 
(Bergvall & Qlofsson [T986l Bergvall et al. 20001 12UU2I Ostlin 
et al., 1999 2001 ). It has unique properties that makes it one of 
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the most relevant local objects for our purpose. The chemical 
abundances are low, the oxygen abundance being ~20% solar 
(relative to the revised solar oxygen abundance by Asplund et 
al. l2()()4b> . 

An estimate of the typical hydrogen mass of a BCG can 
be obtained from the correlation between gas mass and major 
diameter (Gordon & Gottesman 1981}. Haro 11 has a major 
diameter of 11 kpc at /ifl=25 mag arcsec~ 2 indicating that the 
H I mass should be approximately 2xl0 9 Af©. Despite sev- 
eral efforts at VLA, Nancay and Parkes (all unpublished) we 
have not been able to detect H I in this galaxy. We estimate 
the upper limit to M(H I)~ 10 8 M B . The extraordinary strong 
[C II]/ll58/i line, observed with ISO (Bergvall et al.|2000l, 
shows that a large part of the H I gas is located in photodis- 
sociation regions (PDRs). We have estimated the contribution 
from H2 and H II gas to between 10 8 and 10 9 Me each. This 
would leave a remarkably small fraction of the total gas mass in 
atomic form. A possible explanation is that the bulk of the gas 
in the halo that normally is in neutral state has become ionized 
by the starburst. Our Her observations show that the main body 
is indeed surrounded by an extended halo of ionized gas. The 
conditions for LyC photon escape thus appear to be unusually 
favourable. It is interesting to note that the nearby BCG Pox 
186, sometimes called 'ultracompact', is also H I quiet (Begum 
& Chengalur l2005l . 

HST GHRS UV spectroscopic observations of Haro 11 
were carried out by Kunth et al. ( 2003 ) who found a quite com- 
plex absorption system with velocity components of Lya in ab- 
sorption on both sides of the emission line. Thus it appears that 
there are both outflows and infall of gas in the centre direction. 
From the absorption lines Kunth et al. estimate column densi- 
ties and find values as high as log n# ~ 20 ctrT 2 . The aperture 
used is small however (2" circular diameter) and covers only 
the very central part of the galaxy. Our optical spectroscopy 
(Ostlin et al. 2005 ) reveals velocity gradients of both gas and 
stars as large as 200 km s _1 across two of the bright conden- 
sations in the very centre, indicating an extremely high mass 
density in this region. 

2. Observations and reductions 

2.1. IUE observations 

In 1984 we observed the galaxy in low dispersion mode with 
the SWP UV spectrograph of the IUE spacecraft. The aperture 
used was 10"x20". A low resolution spectrum is shown in Fig. 
H Below we will use the slope of the spectrum as an input to 
the derivation of the dust extinction. 

We will also use the IUE spectrum to determine the C/O 
abundance ratio, based on the C III]/il908, and O III]il663 
lines. The strongest line, also representing the dominating ion- 
ization stage is C III] Al 908 while the O III] line is nosier. 
Thus the error in the determination will be significant but it 
still will be useful as a check of the C/O ratio determined from 
the FUSE absorption lines (see sect. 3.1). The atomic collision 
strengths are determined from the electron temperature derived 
by Bergvall & Ostlin (2002) T « 13400 K. From the relative 
line intensities we then obtain N(C + )/N(0 ++ ) ~ 0.3 ± 0.1. The 
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Fig. 1. Spectrum of Haro 1 1 obtained with the IUE space- 
craft in the low resolution SWP mode. The aperture used was 
10"x20". Bad pixels have been flagged with crosses. 



error is dominated by the photon statistics and has been esti- 
mated from the noise in the region around the line. Assuming 
ionization correction factors ICF(C/C + )~ 1.4 (Stasiriska|l990l 
and ICF(0/0 ++ )~ 1.5 (Bergvall & 0stlin l2002l . we then finally 
obtain N(C)/N(O)=0.46!° : [5. 

2.2. FUSE observations 

The data were obtained on Oct 12, 2001, using the FUSE low 
resolution mode LWRS with an aperture of 30"x30". The total 
integration time on target was 16ks. 12ks of the observations 
were carried out during orbital night. The data were reduced 
using CalFUSE v3.0, with default processing parameter for all 
spectral segments. As the background level in FUSE data is 
quite low and hence cannot be reliably measured in detail for 
short observations, the standard FUSE background correction 
is based on a set of template files which are scaled to averages 
over the whole detector in any given observation. These tem- 
plate files consist of coadded data from many, long observa- 
tions of blank sky acquired throughout the mission. The FUSE 
detectors consist of microchannel plates with electronic pixels 
whose exact locations depend on the electrical field structure 
and read-out electronics of the detector (Sahnow et al., 2000). 
Artificial count enhancements occur at the edges of the detec- 
tor due to non-uniformities in the field and have a very steep 
dependence on the distance from the edge. The exact mapping 
of the detector to location in the electronic detector-image de- 
pends on several factors including the total brightness of the de- 
tector illumination. Because of the steepness of this edge effect, 
minor mismatches in the target and background images can 
cause significant errors in the background subtraction (Sahnow, 
2005, private communication). As the SiC spectrum falls very 
close to the lower edge of detector IB, any relative change in 
the strength of this edge effect compared to the rest of the back- 
ground will cause problems in the background correction, par- 
ticularly for weak sources such as Haro 1 1 . To compensate for 
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this effect, we reset the areas and relative weighting used to 
determine the scaling of the background template in the reduc- 
tion of the SiClB data. Without any further modifications the 
derived SiClB spectrum then agrees well with the SiC2A spec- 
trum. The regions of the lowest signal in both the SiClB and 
the LiFl A segments were however significantly below the zero 
level, indicating that the proper background on the whole was 
fainter than the subtracted template background. After extrac- 
tion, calibration and combination of the individual segments, a 
constant was therefore added to the background to compensate 
for this. As we are more strongly constrained by low signal- 
to-noise than spectral resolution, we are not significantly af- 
fected by the slight difference in dispersion solution between 
segments, and the consequent lowering of resolution caused by 
this coaddition. 

3. Results 

Fig. [5] shows the full spectrum with a few of the strongest 
absorption and emission lines indicated. Strong geocoronal 
Lyman emission lines are seen on top of the absorption spec- 
trum. The narrow spike at 1168.7 A is a scattered Hel solar 
emission line. From the absorption lines we can define two 
velocity systems, one originating in the Milky Way and the 
other in Haro 11. Fig. |2] shows the H I 21-cm line spectrum 
in the direction of Haro 11 (1=328°, b=-83°), obtained from 
the Leiden/Dwingeloo survey (Burton & Hartmann 1994 1. Our 
analysis shows that the spectrum in the velocity range -450 
km s _1 ~ 400 km s~'can be approximated with two components 
with small negative velocities relative to the LSR as shown in 
TableU] 
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Fig. 2. The 21-m line in the direction of Haro 1 1 obtained from 
the Leiden/Dwingeloo survey (Burton & Hartmann 1994 1. 
Vlsr is the radial velocity relative to the local standard of rest 
and T A is the antenna temperature. 



The spectral features identified to belong to the Milky Way 
ISM are indicated with hatched vertical lines in Fig. [3] These 
absorption lines have a linewidth of FWHM ~ 30 km s _1 in 



Table 1. Components of the local H I spectrum in the direction 
of Haro 1 1 . The columns show the central velocity relative to 
the local standard of rest, peak antenna temperature, full width 
at half maximum and the derived column density. 
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agreement with the FWHM given in Tabled The narrowness 
of the lines is important since it ensures that we can find the 
proper level of the Lyman continuum of Haro 1 1 between the 
positions of the galactic Lyman absorption lines. The spectrum 
of Haro 1 1 is rich in spectral features and dominated by absorp- 
tion lines from O and B stars. Possible detections of molecular 
hydrogen lines in the Lyman series have also been indicated. 
The heliocentric systemic velocity of Haro 1 1, as defined from 
the optical H II emission (unpublished) and the stellar Call 
triplet absorption lines, weighted from the UV mapping of the 
galaxy (Ostlin et al. 12005b . is v=6180 km s _1 . This velocity 
agrees with the position of the P V Al 118 A (A obs =U4l A) 
line. The P V line originates from hot dwarf stars on the main 
sequence. If these stars are metal-rich and massive they tend 
to show P Cygni profiles while at lower metallicities, as in our 
case, the feature is mainly photospheric and narrow and can be 
used as a systemic velocity fixpoint. Thus we have two inde- 
pendent derivations of the systemic velocity of the stellar pop- 
ulation seen in the far-UV spectral region. Strong deviations 
from this value would indicate in- or outflows of gas. 

The strongest absorption lines, like those of the Lyman se- 
ries, are dark but not saturated as is normally the case for star- 
burst dwarfs (e.g. Heckman et al. 2001 1. This in itself is an indi- 
cation that the optical depth in the LyC is lower than normal for 
this type of galaxies. Many of the absorption lines originating 
in Haro 1 1 are broad and show P Cygni profiles, e.g. the O VI 
/LU032,1038 (redshifted to ,11053,1058 A) lines. Most of these 
lines are probably formed in the photospheres or in the winds of 
O and B stars but there may also be a contribution from more 
global outflows. After a rough correction for the stellar con- 
tribution (see next section), some of the strongest absorption 
features show asymmetries, indicating high-velocity outflows. 
This aspect will be briefly discussed below and more in detail in 
a forthcoming paper, focused on the absorption line spectrum. 

In Fig.0]we display the short wavelength part of the spec- 
trum. This spectral region is covered by two detector segments, 
IB and 2 A, producing two independent spectra SiClB and 
SiC2A. Wavelengths below 917 A are only covered by seg- 
ment IB. This is also the detector segment with the small- 
est effective area of the two. Thus, in this region the signal is 
noisier than at longer wavelengths. What is particularly inter- 
esting from Fig. [4] is that below the Lyman limit of Haro 11 
there clearly seems to be a weak but significant signal in the 
Lyman continuum. Shortwards of the Lyman limit of the local 
rest frame the signal diminishes gradually, as expected from 
the increasing density of the Lyman absorption lines of the lo- 
cal ISM. The emission is seen in the 2-dimensional display of 
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Fig. 3. FUSE spectrum of Haroll. The major absorption features have been indicated with solid lines (Haro 11 optical rest 
frame) and hatched lines (local frame). The spectrum also shows strong geocoronal Lyman emission lines and scattered solar Hel 
at Al 169 A. The data have been rebinned with a boxcar of size 0.26 A. 



the spectral region as a faint narrow noisy structure in both SiC 
channels. We therefore feel confident that the LyC excess is 
not due to calibration problems. The level of the signal how- 
ever is not so well determined, mainly because of uncertainties 
in the calibration, a low signal and possibly influences of auro- 
ral emission lines and emission and absorption lines from the 
Milky Way which we will now discuss. 

Haro 11 is at a galactic latitude of -83 deg. and 
thus interstellar absorption is not a serious problem. From 
Schlegel at al. ( 1998) we obtain A B =0.049 m corresponding 
to Aiy-Hmit ~0.18 m (Mathis 1990 1. The Lyman emission lines 
from the terrestrial airglow are indicated in Fig. [5] For com- 
parison we show in Fig.|S]both the spectrum based on night- 
and-day data and the spectrum based on night data only. We 
can clearly identify the geocoronal Lyman emission lines. 
Shortwards of the Lyman lines, at much lower intensities, C I 
lines may weakly influence the counts (Feldman et al. 200 1}- In 



the following discussion we will use the 900 A rest frame flux 
density of Haro 1 1 as a reference of the LyC level. This choice 
allows comparisons with previous studies of other galaxies 
(e.g. Leitherer et al. 1199 5 Deharveng et al. 120011 Inoue et al. 
2002 ). In the local frame this part of the spectrum has been red- 
shifted to A-920 A. To beat down the noise in the spectrum we 
determine the LyC level over the line-free regions in the inter- 
val 920-925 A, since the LyC signal is not expected to vary 
significantly over such a short interval. After correcting for 
galactic extinction, we obtain f9oo=l.l+0.1xlO" 17 Wm _2 A _1 . 
The uncertainty due to the complexity in the determination of 
the background is probably significantly larger than the formal 
error. In the worst case it could amount to of the order of ~ 50% 
of the signal. 
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Fig. 6. Comparison between night and day data. The lower panel shows the 1 cr noise level, corrected to the binsize used in the 
plot. 




3.1. Spatial distribution of the absorbers. 

We mentioned in the introduction that Haro 1 1 appears to be 
devoid of neutral hydrogen. This could indicate that the LyC 
photons are leaking out from a truncated Stromgren sphere, 
i.e. that the emission is 'density bounded' . The hydrogen ion- 



ization cross section at the Lyman limit is « 8 x 10~ 18 cm 2 . 
This means that the LyC is essentially black at column den- 
sities of a few times 10 17 cirT 2 , if the neutral gas is diffusely 
distributed. As discussed in sect. 1.1 we have derived an upper 
limit of the mass of neutral hydrogen in Haro 11 of 10 8 M & . 
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We can use this number to estimate the maximum column den- 
sity of neutral hydrogen under the assumption that the H I has 
the same spatial distribution as the stars in Haro 1 1 . The stel- 
lar scalelength is hy ~ 2.5 kpc (Bergvall & Ostlin 2002). This 
gives an H I column density in the direction of the centre of 
N(H I)=5xl0 19 cm~ 2 which can be regarded as an upper limit, 
considerably larger than the limiting column density for Lyman 
leakage. In this simple model H I would become transparent at 
4-5 scalelengths (i.e. ~ 10 kpc) from the centre, where the host 
galaxy component starts to dominate. But we know that the 
neutral gas is clumped so we will adapt the simplest alternative 
to the homogeneous model - the 'picket-fence' model. Here it 
is assumed that the LyC radiation is leaking out through holes 
in the ambient gas. In principle we can carry out profile fitting 
to the Lyman lines in absorption to derive the H I column den- 
sities. These lines are however not very useful since there is a 
mixed contribution from OB star photospheres, ISM absorption 
and line emission from the H II regions. We can make an esti- 
mate of the influence of stellar features from models of spec- 
tra of young stellar populations by Robert et al. (2003). These 
models are based on observed stellar spectra over a range of 
metallicities. The metal-poor models are valid for ages below 
15 Myr. We interpolated the model spectra to a metallicity in 
agreement with what we obtained from the optical spectra, i.e. 
an oxygen abundance of ~ 1 8 % solar (based on the solar oxygen 
abundance derived by Asplund et al. 2004b). We also assumed 
a Salpeter initial mass function (IMF) and a constant star for- 
mation rate (SFR). Fig. shows the best fit to these models, 
a burst of the highest age in the library, 14 Myr. The optical 
spectra of the central burst show strong W-R features and, as 
we mention above, P Cygni profiles are seen in the FUSE spec- 
trum. These do not fit to the model spectrum shown in fig. 
A better fit can be obtained if we assume an increasing SFR 
with the UV emission dominated by massive stars of an age of 
a few Myr. Modifications of the standard IMF can further im- 
prove the fit. There may also be a significant contribution to the 
P Cygni component from the supernova generated superwind. 
However, these issues are not of main importance here. 

The influence of young stars in the spectrum is strong. The 
presence of O VI and its P Cygni profile indicates the presence 
of O supergiants (e.g. Gonzalez Delgado et al. 1997 ). Thus the 
burst must be active or stopped not more than 3-4 Myrs ago. 
For a continuous SFR we do not expect the spectrum to change 
much over the expected maximum duration of a typical burst, 
~ 10 8 yr and therefore the comparison in Fig.[7]should be qual- 
itatively correct. The model spectrum has only partly been cor- 
rected for interstellar absorption but the influence on the LyyS 
line is marginal so the line profile should be entirely dominated 
by the stellar photospheric contribution. This is confirmed also 
from a comparison with predicted Ly/3 absorption line strengths 
in models of O star spectra (Gonzalez Delgado et al. 1997). 
Thus it is obvious that these features have a strong influence 
and thus the Lyman lines cannot be used to derive the column 
density without a more careful correction for the influence of 
the stellar features. However, if we can derive column densities 
of low ionization stages of heavy elements with known abun- 
dances relative to hydrogen we can obtain a stricter constraint 
on the H I column density since the atoms and ions in these 
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Fig. 7. The FUSE spectrum of Haro 1 1 (solid line) along with 
the best fitting stellar spectrum (dotted line) from Robert et al. 
(2003). The wavelength scale has been corrected to the rest 
frame of Harol 1, assuming a systemic heliocentric velocity of 
6180 km s _I . The model spectrum is based on a Salpeter IMF, 
a stellar mass range of 1-120 /Vt Q , a metallicity of 18% solar 
and a continuous SFR. The age is 14 Myr. 

stages to a first approximation should coexist with the neutral 
hydrogen. 

We will use a few rather well defined and isolated low ion- 
ization lines, C II ^1036, O I ^1039, Sill /H021, Arl ,11048, 
Fell /U097 and Al 125. We cannot exclude the fact that there 
may be a weak contribution to some of these lines from the 
stellar photospheres but the bulk of the absorption should orig- 
inate in the ISM. The ionization potentials of these lines are 
slightly higher than 1 Ryd. It is known from studies of the local 
interstellar medium that elements with ionization potentials be- 
low that of hydrogen and the ionization potential of the second 
ionization stage above, exist predominantly in the singly ion- 
ized stage where hydrogen is neutral. It is therefore possible 
to use also these lines to obtain good estimates of N(HI) un- 
der certain conditions. The lines of the neutral atoms O and Ar 
can directly be used to derive column densities if the atoms are 
effectively shielded from ionizing photons by the neutral hy- 
drogen. In regions where H I is partially ionized, argon may be 
more fully ionized as a consequence of it high cross section to 
photon absorption. If so, we would underestimate the column 
density based on Ar I. This seems to be confirmed from the 
results in Table|5] Oxygen behaves differently since the ioniza- 
tion fraction is coupled to that of hydrogen via resonant charge- 
exchange reactions (Jenkins et al. 2000). 

The optical depth of a specific species is given by: 

2 

Tie i 

r{X) = jfA 2 N(A) (1) 

m e c 

where /is the oscillator strength. From this we can derive 
the column density in cm" 2 : 

r 3.768 • io 14 r 

N = J N(v)dv = — J T(v)dv (2) 



8 



Nils Bergvall et al.: First detection of Lyman continuum photon escape from a local starburst galaxy 



where A, the wavelength of the line, is given in A and v is 
the velocity in km s _1 . 

We derived the column densities using the Owens software 
obtained from the French FUSE team. The line minimum cor- 
responds to the optically derived radial velocity. The typical 
virial velocity of BCGs of this luminosity is < 100 kms~' . The 
width of the absorption lines are much larger than this. This is 
probably due to outflows from the starburst and to a lesser ex- 
tent also to infall of neutral gas. We have to adopt the fit of the 
line profiles to this condition and the programme we use con- 
veniently allows one to fit lines belonging to several velocity 
components simultaneously. Our approach was to test a num- 
ber of situations in which we assumed the cloud ensemble to 
consist of between 3 and 9 different velocity systems, in all 
cases covering the whole width of the lines. In the cases where 
the number of components were low this was compensated by 
the software by increasing the Doppler broadening parameter, 
b. The b parameter was allowed to vary between 30 and 180 
kms -1 . FUSE has a resolution of ~20 kms -1 and one might 
worry that there could be a significant number of unresolved 
saturated components with low b parameters in the sightline. 
From FUSE curve-of-growth of neutral and low ionization ab- 
sorption line systems in the Magellanic bridge one obtains b 
parameters in the range 15-20 kms -1 (Sembach et al. 1979 
Lehner 2002 1. The H2 lines give smaller b values, indicating a 
stronger confinement and lower turbulence of these clouds. We 
do not expect the conditions in Haro 11 to be more quiet than 
in the Magellanic bridge and therefore we should be safe with a 
lower value of the b parameter of the low ionization lines of 30 
kms -1 . In the example we discuss below, however, we allow the 
parameter to vary freely, utilizing the full resolution of FUSE. 
In one setup of parameters we fixed the relative velocity differ- 
ences between the different components to the same value, so 
that the full width of the lines was covered. We also tested fits 
where both the velocity differences and the line widths were al- 
lowed to vary freely. In all cases we assumed a gas temperature 
of 300K. The absorption lines from all atoms and ions under 
study were included simultaneously in the fits. The optimum 
fits were found to contain between 5-7 components. Ideally 
one would expect the relative strengths of the components to 
be the same for each element. As the components get weaker 
and noiser however, the results become more and more uncer- 
tain. Thus while the central 3-4 components show a good agree- 
ment, the fainter components may deviate significantly from 
one species to the other. The goodness of fit (lowest total x 1 
residual/numbers of degrees of freedom) showed small varia- 
tions between fits to different numbers of subcomponents. As 
an example, one of best results for the carbon line is shown in 

Fig.ia 

From Table |2] we obtain an oxygen to carbon ratio of 
N(C)/N(O)~0.5. This agrees with the value derived from the 
IUE spectra (see sect. 2.1). Thus we feel confident that the pro- 
cedure we applied above gives reliable results. From the fits 
we derived the column densities listed in Table |2 We have 
previously obtained the oxygen abundance of H II regions in 
Haro 11, log(N(0)/N(H))+12=7.9. We will assume that this is 
an upper limit of the oxygen abundance in the neutral region. 
Applying the conversion to hydrogen column densities from 




1056 1058 
X (A) 



1060 



Fig. 8. The best free fit to the observed C II A 1036 line based 
on 5 velocity components with different Doppler widths and 
component separations. The strongest component is centered 
on the systemic velocity as derived from the optical data. 

Table 2. Column densities of low ionization species 



Species 


A 


logN 


0" logN 


O 


1039.2 


15.48 


0.02 


C 


1036.3 


15.19 


0.02 


Si 


1020.7 


15.10 


0.03 


Ar 


1048.2 


13.67 


0.03 


Fe 


1096.9 


14.80 


0.08 




1125.4 


14.59 


0.13 



this number we obtain a lower limit to the hydrogen column 
density, logN(H I)=19.6. We have no corresponding determi- 
nation of the remaining elements but if we assume the rela- 
tive abundances to be the same as the solar (Asplund et al., 
2004a 1 we obtain a lower limit to logN(H I) of between 19.6 
and 20.3 with the lower value for the strong carbon line. These 
relative high values are not consistent with a diffuse distribu- 
tion of H I since it would prohibit leakage. We conclude that 
the photons must be leaking out through transparent windows 
in the gaseous halo in a "picket-fence fashion". This would also 
mean that the estimated column densities are slightly too low, 
due to the dilution. A possible scenario explaining why we ob- 
serve a leakage is that the merging process has allowed young 
star clusters to separate from high density gaseous regions on 
a time scale shorter than the lifetime of the burst, i.e. in a short 
moment of strong gravitational disturbances. 

In the following analysis we will assume that the picket- 
fence model is the best approximation of the situation. 

3.2. Extinction corrections 

In the next section we will derive an estimate of the LyC escape 
fraction based on both the observed Ha flux and the continuum 
fluxes at 900 and 960 A rest wavelength of Haro 11, fgoo and 
fgeo- The fluxes were integrated over spectral regions of widths 
5 and 2 A respectively. As mentioned above, the 900 A contin- 
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uum wavelength was chosen to allow a direct comparison with 
previous studies of Lyman leakage. The 960 A wavelength 
was chosen for spectral evolutionary model comparisons since 
it represents a spectral region reasonably close to the Lyman 
limit but unaffected by Lyman absorption line wings and other 
strong lines. The observed data need to be corrected for galactic 
and internal dust extinction. The galactic extinction in the di- 
rection of Haro 1 1 is A B =0.049 (Schlegel et al. lT998l . The 900 
and 960 A rest wavelengths are redshifted to ~920 and 980 
A in the observed spectrum. At these wavelength Ab=0.049 
corresponds to A 920 =0.177 and A 98 o=0.175 (Mathis [T990t . We 
will now discuss the internal extinction corrections. 

In the picket-fence approximation, as is further detailed in 
the next section, we will assume that the LyC photons are es- 
caping through dust free windows. The only goal for the re- 
maining part of this section is thus to obtain a value of the ex- 
tinction correction to the 960 A region, A960- From the HafHB 
ratio of the central 4"x4", corrected for underlying absorption 
(Bergvall and Ostlin 2002 1, we obtain, assuming the Calzetti 
et al. i 19941 1 extinction law, an internal extinction in B of 
A B =1.08 m +0.07 m . The FUSE aperture is 30"x30", i.e much 
larger than the aperture of the optical spectroscopy. The flux 
inside the 4"x4" aperture is only about 10% of the total flux 
so an estimate of the global extinction based on these data is 
quite uncertain. If we assume that the dust column density is 
proportional to the B band flux we use the small aperture data 
to derive total extinction in B to be A B =0.66'"±0.08 m . The er- 
ror does not include the uncertainty of the model of the dust 
distribution. 

We can also estimate the extinction from the observed slope 
of the UV continuum in Haro 11, using a method developed by 
Calzetti et al. < fl994l [19951 and Meurer et al. ( fl995l) . Calzetti 
et al. ( 1994) investigated 39 starburst and blue compact galax- 
ies to study the interstellar extinction in the UV. The observed 
spectral distribution in this region may be approximated by a 
power law F(/l)oc A 13 . They found that the optical depth obtained 
from the HafHB ratio correlated with the slope B and derived a 
value of B in a dust free environment, j6o=-l. 7. They also found 
that the optical depth calculated from the HafHB ratio was a 
factor of 2 larger than that derived from the continuum under- 
lying the emission lines. This was interpreted as a consequence 
of a systematically higher extinction in young star forming re- 
gions as compared to older regions. 

Later, Meurer et al. ( 1995) found a correlation between B 
and the ratio between the far infrared flux and the UV flux at 
2200 A. This correlation agreed very nicely with a foreground 
dust screen model following the Calzetti et al. extinction law, 
if the galaxies had a similar intrinsic UV continuum distribu- 
tion. The observed slope can therefore be used to derive an esti- 
mate of the extinction by comparing to the "standard" dust free 
slope. The value derived by Meurer et al. Bq--2.5±0.5, is bluer 
than what Calzetti et al. found. This is probably due to the fact 
that the Meurer et al. sample contained more low-luminosity 
starburst dwarfs that were completely dominated by the young 
population while the more massive galaxies in the Calzetti et 
al. sample had a slightly larger contribution from older stars. 
The extinction of Haro 1 1 therefore depends on which value we 
choose. Haro 1 1 is involved in a very intense global starburst 



and is metal poor. On the other hand it is quite luminous and 
most likely has quite a complex star formation history so we 
expect a mixture of stellar populations of different ages. This 
puts Haro 1 1 somewhere in between the Meurer et al. sample 
and the Calzetti et al. sample. We will therefore assume that its 
dust free UV continuum is represented by the weighted mean 
of the Meurer et al. value and the Calzetti et al. value. Both 
have a similar observed scatter around the mean so the value 
we will adopt is Bq--2. 1+0.2. 

The UV slope in Haro 11, obtained from our IUE obser- 
vations, covering the wavelength interval 1250 A-1980 A is 
/?=- 1.4+0.1. (In the picket-fence model it is incorrect to use 
the slope without correcting for the contribution from the 
unobscured part, but this would make a negligible differ- 
ence here.) As an additional check of this value we can use 
the observed flux ratio Ffm/Fuv of Haro 11 and apply the 
model most favoured by Meurer et al. to calculate the pre- 
dicted value of B. F™ is derived from IRAS data and the ap- 
proximation from Lonsdale & Helou (1985). This results in 
F™=1.26-10- 14 (2.58F 60 +Fioo) W m" 2 = 2.74-10" 13 W trT 2 . 
The flux at 2200 A was obtained from an extrapolation of the fit 
to the shorter wavelength region. This gives us Ffi^fFuv-5-9. 
From this value we obtain B-- 1 .6+0. 1 , agreeing nicely with the 
former value we derived. B--1.4 corresponds to an extinction 
in B of Ab=0.83, which is quite close to the value we obtain 
from the HafHB ratio (Ag=0.66). We will assume a value of 
A B =0.7±0.2. 

The final value of the internal extinction was obtained in 
two steps. First we adopted the Calzetti et al. extinction law, rel- 
evant in the wavelength region 1250 A-8000 A , to determine 
the reddening correction at 1250 A. Then we used the recent 
determination of the extinction curve in the SMC by Cartledge 
et al. ( 2005 ) to go from there to 960 A, the region we will use to 
sample the stellar continuum. We thus finally adopt an internal 
extinction A 96 o=2.9 m ±0.3. 

We will now, with the use of spectral evolutionary models, 
derive an estimate of the upper and lower flux escape limits. 

3.3. Estimate of the escape fraction 

3.3.1. Comparison with spectral evolutionary models 

Under the assumption of a simple geometry for the distribution 
of stars, gas and dust, limits on the global LyC escape fraction 
from Haro 11 may be inferred from the observed LyC flux. 
In principle, the LyC leakage could be caused either by holes 
in the gas (e.g. supernova chimneys), by a total gas mass too 
low to form a complete Stromgren sphere, or by a combination 
thereof. Since the high column densities of neutral hydrogen 
derived in Sect. 3.1 disfavour a truncated Stromgren sphere for 
Haro 11, we will here estimate the escape fraction assuming 
that the escape is exclusively due to an isotropic distribution 
of holes, which are assumed to be devoid of gas (neutral and 
ionized) as well as dust. 

Following Deharveng et al. (2001 1, we define the total LyC 
escape fraction / esc as 

, _ £-900,obs 

/esc — ~Z > (J) 

^900,+ 
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where Lgoo.obs represents the observed (escaping) luminosity 
(W A -1 ) at a wavelength corresponding to 900 A in the rest 
system of the galaxy, and L900,* represents the corresponding 
luminosity produced by the stellar component before attenua- 
tion by gas. Under the approximations adopted here, L900,* is 
given by the sum of two components: 



^900,* - ^900,obs + i^OO.absj 



(4) 



where Lgoo.abs represents the luminosity lost through gas ab- 
sorption inside the galaxy. Lgoo.abs is in turn related to the num- 
ber of LyC photons absorbed by gas JVLyCabs through 



i900,abs - A^LyC.abs X 10 



(5) 



where k\ is a model-dependent parameter. Under the assump- 
tion of case B combination, the number of absorbed LyC pho- 
tons can then be related to the extinction-corrected Ha lumi- 
nosity L]i a through 



m/3 



a B 



hc(jua/j mathrmHp) 



-Ha, 



(6) 



where jn a / is the intrinsic Ha/H/3 line ratio, is the Case 

B recombination coefficient and af?„ is the effective recombina- 
nt 

tion coefficient for the Hfi line. These three parameters depend 
on the temperature of the nebula and, to lesser extent, on the 
density. By assuming an electron density of 100 ctrT 3 and fit- 
ting a cubic spline to the ju a / jap, and a^ values given by 
Osterbrock ( 1989), we estimate that for the electron tempera- 
ture of Haro 1 1 (T * 13400 K; Bergvall & Ostlin l2002t : 



^LyC,ab S *8.56xl0 n LH ff . 



(7) 



Combination of the equations above, with the relation L = 
4nD 2 f between luminosity L, Ha flux fya, flux density fx and 
distance D, gives: 



/900,obs 



(8) 



f mhobs +$.56xl0 u - k >f Ha ' 

Constraints on /esc may be derived from by combining the 
estimated l<x error bars on /9oo,obs and fn a with extremum val- 
ues of k\ derived from models. When evaluating the plausible 
range of the k\ parameter, defined as 



N LyC ,abs 

h = log(- ), 

^900,abs 



(9) 



we use the Zackrisson et al. (2001 hereafter Z01) spectral evo- 
lutionary model. With this model, a grid of burst-like spectral 
evolutionary sequences is generated, defined by the model pa- 
rameter values specified in in Table [3] In the case of a stan- 
dard Salpeter IMF (dN/dM oc M~ a \ a = 2.35 with upper mass 
limit M up = 120 M ) the range of k\ values allowed by this 
grid of star formation scenarios becomes k\ = 13.23-13.24. 
When the full range of IMF variations considered (a = 1.35- 
3.35 and M up = 20-120 M ) is used, this constraint relaxes to 
ki = 12.70-13.58. 

When these values of k\ are combined with the previously 
discussed measurements of /9oo,obs = 1 .06+ 0.09 x 10~ 17 W irT 2 
A -1 and / H „ = 4.5 + 0.5 x 10~ 15 W irT 2 , limits on / esc may be 



inferred, resulting in 0.04 < / esc < 0.10 for the standard IMF 
and 0.01 < / esc < 0.11 for the more generous range of IMF 
scenarios. 

These / esc estimates rely heavily on the conversion from 
Lu a to Lgoo.abs, and are only valid under the assumption that the 
Ha line is dominated by photoionization by stars, i.e. that ion- 
ization contribution from shocks or an active galactic nucleus 
are negligible. We have also assumed the dust to be exclusively 
located outside the Stromgren sphere, so that no LyC photons 
are lost due to internal dust extinction. Potential far-UV opac- 
ity sources outside Haro 1 1 , other than dust in the Milky Way, 
have furthermore been neglected. 

A consistency test may be carried out by comparing the ob- 
served slope of the spectrum across the Lyman limit, to the cor- 
responding slope predicted by models. To do this, we measure 
the continuum flux density at a rest wavelength of 960 A, where 
the spectrum should be completely unaffected by gaseous ab- 
sorption. After correction for Galactic extinction, the observed 
flux density at this wavelength is /96o,obs = 7.4 + 0.4 x 10 -17 
W irT 2 A -1 . In sect. 3.2 the correction factor for internal ex- 
tinction was estimated to be yg^o = 14^. Given the assumption 
that the holes in the nebula are completely devoid of both gas 
and dust, and hence that the dust correction applies only to the 
parts of the nebula where there is no LyC leakage, the escape 
fraction may then be estimated from: 



fesc — 



^2/9()0,obs 



y96o/960,obs _ (7960 _ l)fe/900,c 



(10) 



where hi represents the ratio of the stellar component fluxes at 
960 A to 900 A prior to absorption by gas and dust: 

, /960,* 



/900,* 



(ID 



Although the exact value of £2 is not known for Haro 11, a 
lower limit on / esc can be imposed from the minimum &2 pre- 
dicted by the plausible range of star formation scenarios. From 
the model grid defined in Tabled we infer tiling) = 1.78 for 
the standard IMF and tiling) = 1-42 otherwise. This trans- 
lates into / esc > 0.025!° °J 2 and / esc > 0.019!°°^, respectively, 
which is completely consistent with the LyC escape fraction 
derived from the Ha flux. 

To further test the robustness of these results, two other 
spectral evolutionary models - Starburst99 v. 4 (Leitherer et 
al. 119391 and PEGASE.2 (Fioc & Rocca-Volmerange[T333} - 
have been used to assess the model-sensitivity of the k\ and ki 
parameters. We find, that when identical input parameters are 
used in all codes, very similar results are produced, indicating 
that our limits on / esc do not critically depend on the use of any 
specific model. 

4. Discussion 

As we mentioned in the introduction it has been argued that 
star-forming galaxies may yield the required LyC emissivities 
to reionize the universe at z >5-6. This approach has the advan- 
tage over assuming AGNs as sources in that it explains the early 
pollution of heavy elements observed in the IGM at these red- 
shifts (e.g. Cowie et al. 1995 1. Moreover, analysis of the optical 
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Table 3. The grid of Z01 evolutionary sequences. The grid con- 
sists of all possible combinations of the parameter values listed 
below. Each evolutionary sequence runs from ages of to 100 
Myr in short time steps. 

IMF, a 1.35, 2.35,3.35 
M up 20, 70, 120 M Q 
SFH elO, clOO 

0.001 

IMF: dN/dM cc Mr" 

M up : The upper mass limit of the IMF 

SFH=Star formation history. 

c=Constant star formation rate during the subsequent 
number of Myr. 

e=Exponentially declining star formation rate, 
SFR(?) oc exp(-f/r), with an e-folding decay rate (t) 
equal to the subsequent number of Myr. 
Z t : The metallicity of the stellar population 



depth at the H I and Hell Lya continuum breaks in different en- 
vironments at z <4 indicates that the major ionization source in 
dense regions is QSOs while in void regions the radiation ap- 
pears to be softer, giving room for a substantial contribution 
from starburst galaxies. 

Several semi-analytical methods have been used to model 
the reionization epoch. Due to the strong inhomogeneities of 
the radiation field in the early stages of the reionization, the 
computations are however quite demanding and complex (see 
Razoumov et al. 2002 for a summary of the technical prob- 
lems) and realistic simulations have not been possible until 
recently. In 2003, Sokasian et al. (2003) used a fast code ca- 
pable of exploring the impact of low-mass starburst galaxies 
in high spatial resolution. They focused on the evolution of 
the interaction between the IGM and starforming galaxies of 
normal properties, e.g. similar to Haro 11, in the redshift in- 
terval z ~20 to z ~6. One of the results was that the previ- 
ous estimates of the required photon flux rate needed to com- 
plete the reionization was overestimated by 63% if assumed 
that only galaxies with masses > 10 9 h _1 A1o would account 
for the reionization. By including sources with masses down 
to ~4xl0 7 Ir'yVlo the impact of the large population of dwarf 
starburst galaxies could be realistically assessed. The calcula- 
tions demonstrated that the soft ionization sources of type star- 
burst dwarfs give the required photon flux provided that the es- 
cape fraction is / esc >0.20. This number is twice the upper limit 
of /esc.gas that we derived for Haro 1 1 . The estimate is how- 
ever based on a normal stellar population. At redshifts z>10 
pop III stars could dominate the population (e.g. Choudhury & 
Ferrara 2005 ). In these stars the emissivity in the UV might be 
larger than in pop II by a factor of 2 (Tumlinson & Shull 2000 
Bromm et al. 2001 ). Thus, although / esc would not change, the 
number of produced photons per unit mass would increase sig- 
nificantly. Thus it seems feasible that, at least at z>10, dwarf 
starbursts of the mass of Haro 1 1 and lower could provide a 
substantial contribution to the early phase of the reionization 
process. 



5. Conclusions 

We report on the discovery of a weak signal in the Lyman con- 
tinuum of the luminous starburst galaxy Haro 1 1 . It is the first 
time Lyman continuum leakage is found in a galaxy in the local 
environment. Using the derived column densities of low ion- 
ization species we estimate the column density of neutral hy- 
drogen to be ~ 10 19 cm~ 2 , which is too high to allow Lyman 
continuum photon escape if the gas is diffusely distributed. 
We conclude that the Lyman continuum radiation is escaping 
through transparent windows of the ISM. Assuming a Salpeter 
IMF, we estimate a total escape fraction of 0.04< / esc <0.10. 
More extreme IMFs allow a lower escape fraction but the major 
uncertainty is due to the problematic background subtraction 
at short wavelengths. The upper limit is below the model es- 
timates of the minimum escape rate required for reionization. 
Still, this single case demonstrates that a significant fraction 
of the required Lyman continuum photons may be produced 
by starburst dwarfs, in particular if we also consider that Pop 
III stars, that should dominate the first generation galaxies, are 
more efficient than Pop II stars in producing Lyman continuum 
photons. 
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